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Abstract. Let f! be a bounded domain in R'^, N > 2, with smooth boundary dU. We 
construct positive weak solutions of the problem Au + uP = in Q, which vanish in suitable 
trace sense on dfl, but which are singular at prescribed single points if p is equal or slightly 
above \ ■ Similar constructions are carried out for solutions which are singular on any given 
embedded submanifold of dfl of dimension 0<fc<7V — 2, ifp equals or it is slightly above 
and even on countable families of these objects, dense on a given closed set. The 
role of this exponent, first discovered by Brezis and Turner T for boundary regularity when 
p < , parallels that of p = for interior singularities. 



1. Introduction and statement of main results 

Let ri be a bounded domain in M^, with smooth boundary dil. A model of nonlinear elliptic 
boundary value problem is the classical Lane-Emden-Fowler equation, 

' Au + uP = in n 

u > in fl (1.1) 

u = on d^l 

where p > 1. We are interested in finding solutions to this problem which are smooth in and 
equal to almost everywhere on with respect to surface measure. More precisely, we want 
to study solutions to problem (|l.ip that satisfy the boundary condition in a suitable trace sense, 
while not necessarily in a continuous fashion. 

Following Brezis & Turner [T] and Quittner & Souplet [llj . we say that a positive function 
u G C°°(r2) is a very weak solution of problem (jl.ip if 



and 



u, uPdist {x,dn) e L\n) 



{u Av + uPv)dx = Q for all v G C'^{n) with v = on dQ. 



From the results in [TJ [TT] , it follows that if p satisfies the constraint 

N + 1 



1 <p < 



N - 1 



(1.2) 
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then a very weak solution u is actually in i/g (fJ), and it is a weak solution in the usual variational 
sense: 

u e H^{n), I (Vu Vw -uPv)dx^O for all v G H^{n). 
Jn 

Elliptic regularity then yields u £ C^{0.), so that u solves classically. As it is well-known, a 
constrained minimization procedure involving Sobolev's embedding implies existence of a weak- 
variational solution to (jl.ip for 1 < p < A natural question is then whether very weak 

solutions of (jl.ip are classical within a broader range of exponents than (|1.2p . Partially answering 
this question negatively, Souplet (12j constructed an example of a positive function a G 
such that Problem (jl.ip . with u'' replaced by a{x)u'P for p > has a very weak solution which 

is unbounded^ developing a point singularity on the boundary. 

The exponent p = is thus critical in what concerns to boundary regularity for very 

weak solutions. The aim of this paper is to construct solutions to Problem (jl.ip with prescribed 
singularities on the boundary. To state an important special case of our main results we need a 
definition: 

Deflnition 1.1. Let u{x) be a function defined in D, and xq G dQ. We say that 

u{x) £ as X Xq non-tang entially 

7T 

lim u{x) — £ for all a G [0, — ), 

Tq (a;o)3a;— >a;o 2 

where Ta{xo) denotes the cone with vertex ^i, and angle a with respect to its axis, the inner 
normal to dfl at xq. 

We have the validity of the following result. 

Theorem 1.1. There exists a number p[^ > such that if p satisfies 

N+1 

J^<P<PN, 

then the following holds: given points C21 ■ ■ • j Cfc ^ 9^1, there exists a very weak solution u to 
problem (jJ.jp such that m G C^(f2 \ {^1, ■ ■ ■ ,^k}) o.'^d 

u{x) —>■ +00 as X non-tang entially, for all i ~ I, . . . , k. 

The study of the behavior near an isolated boundary singularity of any positive solution of 
p.ip when when the exponent p > was recently achieved by Bidaut-Veron-Ponce-Vcron in 
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1.1. The parallel with p = and interior singularities. The role of the exponent p = 
parallels that oi p = j^z^ for solutions to problem p.ip with interior singularities. Let us 
recall that if u e LP{fl) is a positive distributional solution of (|l.ip and 1 < p < j^z^, then u is 
smooth in fi. On the other hand, for p > -^^^j' distributional solutions of (jl.ip with prescribed 
interior singularities are built in [71 [9l [TOl [H [8l [5l [6j . Basic cells in those constructions are radially 
symmetric singular solutions u = u{\x\) for the equation 

Au + uP^O. (1.3) 

Whenever p > -J^g i t^^*^ function 



Uq{\x\) = Cp^N \x\ , Cp^N 



(TV -2 -) 



p-1 p 



(1.4) 



is a explicit singular solution of (|1.3p in \ {0}. If, in addition, jfz2 < P < "w^i phase plane 
analysis for the ODE corresponding to radial solutions of (|1.3p . yields existence of a singular 
positive solution ui which connects the behavior of uq near the origin with fast decay at infinity, 

ui{\x\) ^Cp^n\x\~^{1 + o{1)) as a; ^ 0, (1.5) 

uii\x\) ^ \x\-'^^-^\l + 0(1)) as|a;|^+oo, (1.6) 

(note that N — 2 > j^)- The scalings ux{r) — Xp-^ ui{Xr) with A > are then solutions of 
(jl.3p that have the same behavior near the origin but which become very small as A ^ 0"*" on 
any compact subset of \ {0}. Thus, given points 

the function 

fc 

u,(x) = ^ux{\x - ^j|) 

i=l 

constitutes a "good approximation" for small A > to a singular solution of Problem (|l.ip . 
Linear theory and perturbation arguments lead to establish the presence of an actual solution to 
(jl.ip near -u,, see [8]. When p = a similar construction can be carried out, see [10]. Basic 
cell ui corresponds in this case to a positive radial solution ui of equation (|1.3p in -6(0, 1) with 

uii\x\) ^ CN\x\~''^~^Uog{l/\x\)-^{l + o{l)) asx^O. (1.7) 

N — 2 

In this case the scalings ux{x) — A^~ui(Ax) have the same behavior as ui at the origin, and 
they approach zero as A 0+, uniformly on compact subsets of \ {0}. 

1.2. The basic cells: singular solutions on a half-space. In the construction of the solutions 
predicted by Theorem 11.11 we will follow a scheme similar to that described above for interior 
singularities. Basic cells will now be positive solutions of equation (|1.3p defined on the half-space, 

R^ := {x={xi,...,xn) / xn>0} 
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which vanish on its boundary, with a singularity at the origin. Such solutions are of course not 
radial, and ODE analysis does not apply. Thus, we consider the following two problems: 



Au - 



> 



\{0} 



(1.8) 



\{0}, 



for p > , and 



At 



yJV-1 



U > 



(1.9) 



on 9RfJ! nB+\{0}, 



where B4 



pN 



n 5(0,1). 



Our purpose is to find families of solutions ux with analogous behavior to the radial singular 
ones previously described. Let us consider first the case p > The role of the explicit radial 

solution uo in (II. 4|) is now played by one found by separation of variables: Let us denote by S^~^ 
the half sphere 

(01,...,0Ar)e^^-V^W>O}. 

of the form 



5^-1 := {9 

Looking for a solution of problem ((T 

uo{x) = r 

we arrive at the problem on the half sphere, 

(A5«-i + iV - 1) 0p - f±i (TV - f±i 



4>p{(^n), r = \x\, = -r—r 



> 



(1.10) 



(1.11) 



Since — 1 is the first eigenvalue 



Here A5N-1 designates the Laplace-Beltrami operator in 
of — A51V-1 under Dirichlet boundary conditions, with eigenfunction On, in the considered range 
N — > 0, an application of the mountain pass lemma yields existence of a solution to this 
problem, provided that, additionally, p is subcritical in dimension A^ — 1, namely p < ^35. When 
p tends from above to this solution ceases to exist by uniform vanishing. Alternatively, in 

this regime, a standard application of Crandall-Rabinowitz local bifurcation theorem yields that 
this solution defines a continuous branch in p with asymptotic behavior 

M^n) = cn{N-^)^9n{1 + o{1)), aspi^^. (1.12) 

Nevertheless, the function uq does not suffice for the construction of approximate profiles for those 
of Theorem 1 1 . 1 1 since it is "too large" at infinity. We need an analogue of the radial function ui in 
p.5p - (|1.6p . namely one that behaves like uq near the origin but having fast decay. A "connection" 
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between uo with Poisson's kernel xn/Ix]'^ does indeed exist provided that p is sufficiently close 
to as the following result states. 

Proposition 1.1. There exists a number pfq > ^~^\ , such that for all jfZ^ < p < p^, there 
exists a solution ui{x) to problem ([773) such that 

2 

ui{x) = \x\'~ (j)p{xN/\x\){l + o{l)) as x ^ 0, 
where 4>p solves and 

ui{x) = \x\^^ xiq (1 + o(l)) as \x\ +00 . 

This solution has indeed "fast decay" since iV — 1 > Observe then that the scalings 

u\{x) = Ap^ui(Ax) define a family of solutions to Problem (jl.Sp which have a common, A- 
independent behavior at the origin, but which vanishes uniformly as A —> 0, on compact subsets 
of \{0}. 

When p = jfzr^ there is no solution to problem (|l.lip and thus separation of variables fails. 
On the other hand, we have an exact analogue of the radial solutions ui in (|1.7p . as described by 
the following result. 

Proposition 1.2. There exists a solution ui of Problem \1.9^ such that 

ui{x) = CN \x\^^ log(l/|a;|)^^ xn (1 + o(l)) as x ^ 0. 

We observe that in this case the functions u\{x) = \^~^ui{\x) satisfy that u\{x) uni- 
formly as A ^ 0+ on compact subsets of M.^ \ {0}. 

1.3. Solutions with prescribed singular set: general statements. In reality, the profiles 
given by the above results can also be used to approximate solutions to Problem (|l.ip whose 
singular set is a fc dimensional submanifold of 9f2 with 1 < fc < — 2. For instance, if wi(x'), 
x' G M.^^^ is the solution of (jl.Sp given by Proposition II . II for p close from above to ^E^rTi then 
u{x) = ui{x') solves the same problem in M.^, now with singular set given by a fc-dimensional 
subspace. This is the content of the following result, more general than Theorem II. 1[ whose 
analogue for interior singularities was found in |10|, [8] . 



Theorem 1.2. Let < k < N ^ 2 and let PN-k be the number given by Provosition \l.l\ with N 
replaced by N — k. Given p with 

N-k + 1 

and a k-dimensional submanifold S embedded in dfl, there exist infinitely many (very) weak 
solutions to problem jp such that u g C^(f2 \ S'), and 

u{x) +00 as X Xq non-tang entially, for all xq G S. 
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When fc = 0, we agree that 5 is a finite set of isolated points, so that Theorem II. H is recovered. 
In reahty, the solutions found arise as continua, depending on as many real parameters as number 
of points lie in S. When A; > 1, the solutions we construct are infinite dimensional families. The 
construction actually allows much more: For instance, when p = n-i ' number of points of 
the singular set can be taken to infinity, to total a dense subset of any given closed set A of Oil, 
which can be properly called its singular set. In fact, since the solutions we are interested in 
are smooth in fl, it is natural to define the singular set of a very weak solution u of (jl.ip as the 
complement in of the set of points x G d^l in a neighborhood of which u is smooth. Observe 
that, by definition, the singular set of u is a closed subset of dil. We have the validity of the 
following general result. 

Theorem 1.3. Let < k < N — 2, and ^Zk-i ^ P < PN-k- Let us consider a nonempty closed 
subset A of dVL, which contains a sequence of k- dimensional embedded submanifolds Si, i G N, 
which are also disjoint and satisfy that S :— UiSi is dense in A. Then, there exists a positive 
very weak solution of Problem <\1.1^ whose singular set is exactly A, and such that 

u{x) — s- +0O as X ^ Xq non-tang entially, for all xq £ S, 

and 

u{x) — > as X — > Xo non-tang entially, for all Xq £ dU, \ S. 

This last result and the underlying construction have interesting consequences: for instance, 
for p larger than but close enough to , there are infinitely many very weak solutions of (jl.ip 

whose singular set is any prescribed closed subset of d^, but such that u S Wq''^{Vl) for any 
1 < g < iV f^TT. Therefore, even though u is not identically equal to at each point of dil^ we can 
say that u = on dVl in appropriate sense of traces. 

The proof of these results relies on two basic ingredients: one is the construction of the basic 
cells of Propositions II. Il and fL2l which we carry our in §2. The other ingredient is the analysis 
of invertibility of Laplace's operator for right hand sides that involve singular behavior near a 
point or an embedded manifold of the boundary. After this analysis, which is carried out in §3, 
the proof of Theorem 1 1 . 21 then follows from a fixed point argument. The result of Theorem ll.3l is 
a consequence of an inductive construction taken to the limit under suitable control. 

2. The half-space case: proofs of Propositions 11.11 and 11.21 
It is natural and convenient to look for solutions of (jl.Sp or (|1.9p of the form 

_ 2 

u{x) — \x\ T-'^ 4>{~\og\x\,x/\x\), 
so that the equation Au + — Q reads in terms of 4>{t, 9) , t £ R, 6 £ S'^'^, as 

9,^0 -(iV- 2^) dt<t>~^ {n-'^) cl^+{As«-.+N-l)cj^ + r = 0. (2.1) 
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2.1. Proof of Proposition 11.21 When p = in the language ol equation (|2.ip . problem 
(|1.9p becomes 

d?(f> + Ndt<l>+{AsN-i+N-l) (l) + <j)^ = in {t,,oo) x S^''^ 

> in (<,,oo)x5^-i (2.2) 

= on (t,,oo) X 95+^^ 

We allow here > to be a parameter, which we will choose later to be large. To get a solution 
of problem p.9p we actually need = 0, but this is simply achieved by a translation of cj) in the 
variable. 

The idea is now to look for a solution of this equation of the form 

(bit, 9) = an t-'" ipiie) + d), (2.3) 

where aN,bN are positive constants to be fixed below, and ipi denotes the eigenfunction of 
— Agjv-i associated to the eigenvalue — 1 and normalized so that its i^-norm is equal to 
1. Explicitly, 



(/. 



N-i 9\j da 



1/2- 



'N ' 

\ ■ -t- 

Whcn substituting the function cj) — ajsrt^^'^ as an approximation for a solution of equation 

(|2.2p . we see that for large t, the main order term in the error created is the function 

E{t,0) := NaNbNt-''-'^ ipiie) ~ |ajv t"^" (^i(6')|^. 
We make the following choice for the numbers un and Bn- 

2N 

ifil'-' da 



On — o,N 



N{N - 1) y5«-i 

This election achieves the L^-orthogonality of E to ipi for all t, namely 

/ E(t,-)ipida ^0 foraUt>i*. 

We fix these values in what follows. In terms of i/j in (|2.3p . equation l|2.2p now reads 

d^^p + N dti^ + {AsN-i + N - 1) ij + bN (&Ar + 1) aw 

NbNaNt'^'' -^i + \aNt^ " = 0, 

^- = on (t*,oo) X aS'+"^ 

We further decompose 

^{t,6)^ h{t)^^{9)+Mt,6) (2.4) 
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where tpi{t,6) satisfies 



ipi {t, •) ifidcr = Q for all t > t^,. 



The equation we have to solve then reduces to the coupled system in /2) given by 

({df + Ndt + iAsN-i+N-l))^bi = A^i(V'i,/2) 



(2.5) 



^1 = on {t^,oo) X dS_^_ 



where 



iVl(V'l,/2 



^2(^1, /2 



N(N-l) 

— ^1 — -aN fi 



( 



N + l N + 1 



r 



\aNt 2 ip^+Tp^ + f^ip^\N-i ^\a^t- 

N + 3 



aM t 



( \aNt' 



(2.6) 



</5i + i/'i + /2<Pi I - law i 

iV + l 



1-JV , W + 1 



a""' i 
Ar_i "at f 1 t . 



r/2)¥'irfCT. 

Here n-*^ denotes the L^^-orthogonal projection over the orthogonal complement to (^i, namely 



h(t,-) ipi da, 



n^{h) = h{t,9) - ipi{0) 



and ip is given by (12. 4p . The logic in the resolution of problem (|2.5p is simple: we look for a 

1 — JV 

solution ("01, /2) which is small compared with (pi. We will construct inverses to the linear 
operators defined by the left hand sides of the equations in ()2.5p with suitable bounds that allow, 
for sufhciently large i*, the resolution of the system via contraction mapping principle. Observe 
that so far we have not imposed boundary conditions at t = . We will invert the linear operator 
in tpi, for right hand sides i^-orthogonal to ipi for all t, imposing Dirichlet boundary condition 
at t = . The choice of inverse for the ODE operator in /2 will be basically explicit, and will not 
require imposing boundary conditions. The natural environment to carry out these inversions is 
L°°-weighted spaces. In the next two lemmas we construct these inverses. Thus we consider the 
linear problems 

'{df + Ndt + {AsN-i+N-l))^ = h in {U,oo) x S^'\ 



on diiU.oo) X S^-^) 



for h such that 



JgN-i 'ip{t, ■) ipi da = for alH > i*, 



h{t, •) ipida = for all t > 



(2.7) 



(2.8) 
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and 



[d^ + Nd, + ^^^^) f^g in (t,,^). (2.9) 



We have the vahdity of the foUowing results. 

Lemma 2.1. There exists a constant c > such that the following holds: Given cr > 0, there 
is a t„ > 0, with t^ > if a > and t„ — Q if a — Q, such that, for all t^ > t„, and all 
h e C''((t*, cx)) X S^~^) that satisfies (|g.<gp and f^h S L°°((t,, oo) X S^~^), there exists a solution 
ip — Ti{h) of problem ^2.7^ , which defines a linear operator of h and satisfies the estimate 

Lemma 2.2. Given a > -^^-^7 there exist numbers t„^ c„ > such that for all > t„ and all 
g G C'^((t*, oo)) satisfying t" g G oo)), there exists a solution f = T2{g) of equation ^2.9^ . 

which defines a linear operator of g and satisfies the estimate 

Before proceeding into the proofs of these lemmas, let us conclude the result. 

Conclusion of the proof of Proposition 11.21 Let us fix in the above lemmas any number a 
such that 

N-l ^ „ ^ jV+1 

and > to-. We obtain a solution of problem ()2.5|) if (t/'i, fi) solves the fixed point problem 

(^l,/2) - X(7/^l,/2):=(Ti(7Vi(7^i,/2)),r2(7V2(V'l,/2))), (2.f0) 

in the space of functions 

(V-,/) G C°([t„ oo) X 5^-1) X C°([t„oo)) 

for which the norm 

||(^,/)|U= lirVlloo+AipVlloo 

is finite. Here < f is a positive number which we will fix later. Ti, T2 are the operators 
predicted by Lemmas 12.11 and 12.21 It is directly checked that we have the pointwise estimates 

|iVi(^l,/2)| < A[t-^+t-l|^l|+<-l|/2|], 

(2.11) 

|A^2(^l,/2)| < A[t-^+t-l|^i|+t^|/2|2 + |/2|^], 

where A depends only on A^, whenever IKV'i/)!^ ^ It follows that, 

P"iVl(^l,/2)||oo < A\t, ^ +C'l|i"^l||oo+C'l|i"/2||oo], 

(2.12) 

iV + l N -\ „ 1 2o- 

pl+"iV2(Vl,/2)||oo < A\tl - +P^Vl||oo + (t*^ +t* "-')P"/2||oo]. 

These estimates, together with Lemmas 12.11 and 12.21 yield that if is chosen sufHciently small, 
depending only on a and A^, and t* is taken sufficiently large, then the operator M. applies the 
ball |l(i/'j/)IU ^ M into itself. A similar estimates shows that, also, is a contraction mapping 
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with this norm inside this region. Hence there is a fixed point (ipi, /2) in this baU. The solution 
obtained this way renders the function 

0) ^aNt-^ ipiie) + f2{t) iplie) + Mt, 0) 

positive in (f*, +oo) x S^^^, and it is then a solution of problem (|2.2p . This completes the proof 
of Proposition [L^ □ 

Next we carry out the proofs of the lemmas. 

Proof of Lemma 12.11 . Let us consider first the case u = 0, so that h is bounded. With 
no loss of generality, we also assume = 0. We see then that problem (|2.7|) has at most one 
bounded solution. This can be shown for instance expanding a bounded solution of the equation 
with ft, = in eigenfunctions of the Laplace-Beltrami operator with zero boundary conditions 
on S^^^ . The coefficients in this expansion will be functions of t which correspond to bounded 
solution of certain homogeneous ODE's which only have the zero solution. Thus, we only have 
to prove existence. To do so, let us consider, for any given number t2 > 0, the problem 

idf + Ndt + iAsN-i+N-l))^j = h in (0,t2) X 5^^"! 

(2.13) 

V' = on a( (0,t2) X 5^-1). 

This problem is uniquely solvable since it is just a rephrasing of a Dirichlet problem for the 
Laplacian in a half-annular region. Let us denote by "0 = "0*2 its unique solution. Since, by 
assumption, h{t, ■) is L^-orthogonal to ipi for all t € (0,t2), so is ip. 

It suffices to check that there exists a constant c > independent of i2 > 1 such that 

ll'/'llL~([0,t2]xS^-i) ^ c|l^llL-([0,t2]xS^-i)- (2-14) 

Indeed, assuming this estimate is already proven, we use elliptic estimates together with Ascoli's 
theorem to show that, as t2 tends to oo, the sequence of functions ipt2 converges uniformly to a 
function solution of (|2.7p which satisfies 

Elliptic estimates then imply that 

ll^^llL-([0,oo)xSf-i) + IIV^IIl~([0,oo)xS^-1) ^ Co II /^ll L~ ([0,oo) x S^-^) • (2-15) 

The orthogonality conditions on -0 pass certainly to the limit, and existence of a solution with 
the desired properties thus follows. It remains to prove the uniform estimate (|2.14p . We argue 
by contradiction. Since the result is certainly true when t2 remains bounded, we assume that 
there exists a sequence t2 = ^2.1 tending to 00, functions h — hi and 4'i corresponding solutions 
to problem (|2.13p for which 

ll'^»lli~([0,t2,,]xSf-i) = 1 '^lld .1™^ ll'**llL~([0,t2,,]xSf -^1) = 0- 

We choose U G (0,^2,0 where HV'illiooQj^ . ]xs""^) achieved and define 
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Using elliptic estimates together with Ascoli's theorem, we can extract from {'ipi)i some subse- 
quence which converges uniformly on compact sets to ip, a bounded solution of 

[d^ + Ndt + {AsN-i+N -l))i>^0 (2.16) 

which is either defined on [0, oo) xS";^^^, on (— cx), 0] x 5;^^^ or on (— cxd, cxd) x S*;^^^. Furthermore, 

II^I|l~=1 (2.17) 

with having boundary data. Furthermore ipit, ■) is L^-orthogonal to (pi, for all t. Eigenfunc- 
tion decomposition of ^{t, ■) for the Laplace-Beltrami operator yields that there is non nontrivial 
bounded solution of (|2.16p and this contradicts (|2.17p . This completes the proof of the uniform 
estimate, and thus existence of a unique bounded solution of p.Tp with the desired estimate 
follows. This solution of course defines a linear operator on bounded h. 

To establish the result for a > and > sufficiently large, let us write 

h^r^h and ^ t^'^ 

so that h is bounded. (12. 7p reduces to 

(a2 + 7vat + (As«-i + 7v-i))v3+(^^-^) i,-^-^dti> = h (2.18) 

We can estimate 

II f cr{cr+l) Ncr\ 7 2a ^ J,, , 
II \— r)^-- ^*V^llL-([t.,oo)xS~-i) < 

M (l|V^||ioo([t^,o^)xs«-i) + \\i>\\L^lXu,c^)y.Sl--')) 

where [i can be taken as small as we wish, after choosing > t„ with t„ large enough. The 
resolution of (|2.18p with the desired bound then follows from that of (|2.7p with cr = together 
with a direct linear perturbation argument. This finishes the proof. □ 

Proof of Lemma 12.21 . Observe that a right inverse for the operator &l + N dt on [i*,oo) is 
given by 

/oo fC, 
e-^^y^ e''^g{s)dsdt 

One checks that 

G{g)\\L^((t,, +00)) < ]vj^ (1 - ^1^) P^^'^5llL~((t.,+oo)) 
provided Nt^ — 1 — cr > 0. This follows at once from the computation 

J*. Jt, 
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and hence 

1 cr+l \ / Ns -cr-l J ^ 1 Nt,-a-l 

Jt ^■^^jv^ t 

The result of the lemma then follows from a simple linear perturbation argument, provided t, > 
is chosen so that 

2 II* ^ Lr(5j||L~((U,+oo)) < 2 II* 5l|L~((t.,+oo)), 

and the result is concluded. □ 

2.2. Proof of Proposition [TTTl Recall that, when p G (w-tt, the Mountain Pass Lemma 
yields the existence oi cpp, a nontrivial positive solution of (jl.lip . This solution then induces a 
solution 

_ 2 

Up{x) := \x\ f-i (j)p{x/\x\), 

of Problem ()1.8|) . for which this time we emphasize its dependence on p. We have to show that 
there exists a solution of (jl.Sp which is asymptotic to Up near and it is asymptotic to 

Uoa(X) ■■= \x\ Xn 

at infinity. Let us consider a smooth cut-off function x which is equal to in -Bi(O) and it is 
identically equal to 1 in R^\_B2(0). We will consider the function %) as a first approximation 

for the solution we are looking for. Since, we recall, (f)p approaches zero uniformly as p [ , 
then the same is true for uo.p away from the origin. The result of the proposition relies on a 
perturbation procedure, and this is the reason why we can only show the for exponents p close 
to tnt^t. To carry out this scheme, we shall build a right inverse for the Laplacian relative to the 
following doubly weighted space: 



Definition 2.1. Given 6,5' e M, the space L'^s,{R^) is defined to he the space of functions 

' oo / 



u G L^^(M^ ) for which the following norm 



ir"iiL~,(R^') = 111^' 

is finite 



X\ " M||L~(i3+(l)) + |||a;| w|lL=o(Riv_B^(i)) 



Hence 5 controls the behavior of the function near and 5' the behavior of the function near 
infinity. Let us consider the problem 

Au ^ |a;|-2/ inM^ 

(2.19) 

u = on \ {0}. 

We have the validity of the following result. 

Lemma 2.3. Let (5 £ (1 — iV, 1) and 5' G {~N, 1 — N) he given. There is a constant c > such 
that for each f G L^g,{M.^), there exists a solution u — G{f) of problem i\2.19^ . which defines a 
linear operator in f and can be decomposed as 

u = u + axuoo 
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with 

|a| + ||u||L~, <c||/lk~,. 



Proof. Let us observe that S {N — 2 + S) < N — 1 precisely when 5 ^ (1 — N,l). Therefore, we 
can define (/s* = ipN.s to be the unique, positive solution of 

-{AsN-i+d{S + N^2))(p, = 1 mS^-^ 

(p, = ondS^-\ 

A direct computation shows that 

- \x\''AKN{\xfip,{9))^c\xf. (2.20) 
Assume that / e Lg^g,{M.^). Given ri < 1 < r2, we can solve the equation jxpAu = / in 
n {Br2 — Bri), with boundary conditions. We use the function x i — > f*{x) \x\^ as a barrier 
to prove the pointwise estimate 

\U\ < c||/||ioo_^,(RiV) \xf 

in n {Br2 — -Sri)- Furthermore, given (5 e (1 — TV, 1) we can use the function x i — > f*{x) \x\^ 
as a barrier to prove the estimate 

\u\ < c||/||i»,(Mi^) \xf 

in M^n(B^2 ~Br,). 

We use elliptic regularity theory as well as Ascoli's theorem to pass to the limit as ri tends 
to and r2 tends to cxo. We obtain a solution u of problem (|2.19p which satisfies the pointwise 
estimates 

l"l <C||/|L»,(K«)|X|* 

in (M^ n Bi) - {0}, and 

\U\ < c|l/||ioc^^(RiV) \xf 

in n (M^ — Bi). Finally, the decomposition of the solution u at infinity into 

u = u + a Uoo 

where the function u satisfies 

1^1 <c||/|L~,(R«)kl'' 

follows easily from Green's representation formula. Moreover, we can directly compute the value 
of a. Indeed, integration of the equation over _B+ := M.^ n S,- yields, for r large enough, 

f\x\-^dx^ I druda^ I {dru){r9)r^-^da-a{N -I) [ On da 
dB+ Js1~^ Js"-' 



+ 



Passing to the limit as r tends to c», we obtain the identity 



a{N-l) eNda = - f\x\-^dx, (2.21) 
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and the proof is concluded. □ 

Conclusion of the proof of Proposition [T7ll To find a solution of problem (jl.Sp . we write 

w = (1 - x)up + V, 

where % is a smooth cut-off function which is equal to in Bi and identically equal to in 
- i?2- Let us fix numbers S e {1 - N,l), S' G {-N,l- N) and let G be the operator defined 
in Lemma l2.3l Then, we obtain a solution with the required properties if v solves the fixed point 
problem 

V - -G (|a;p(A(l ~ x) Up + li^ - x) Up + vH) (2-22) 
in the space Lf^gidRf - {0}) © Spanjxuoo}), and (1 - x)"p + v > 0. 
Let us observe that there exists a constant cq — c{N, 5, S') > such that 

IlkP (A((l - X) Up) + ((1 - x) %r)llL~ ,(K«-{o}) < Co Il'^pllc2(sf-i)- (2.23) 
Indeed, since Aup + = we have 

A((I - x) Sp) + ((1 - X) upf = - Ax Up - 2Vx • Vup + (1 - X - (1 - xY) 
and the estimate follows at once. 

On the other hand, if we assume that p is sufficiently close to from above, we have that 

||^p||c2(sJV-l) < 1, 

and also 

6>-^ and 5' > p{l- N) + 2. 

Under these constraints, it is not hard to check the existence of a constant c — c{N, S, 6') > 
such 1jh.3j1j 

|||X|2(|(I -x)Up + V2\P - 1(1 -X)up + Vi\P)\\L^ru--{0}) 

(2.24) 

< c||(/)p||c2(s«-i) 11^2 - wi|lioo_^,(Riv_{o})eSpan{x«oo} 
for aU V2,vi e i;f^-,(R^ - {0}) Span{xMoo} satisfying 

,((R5'-{o})eSpan{xu„}) - 2 co ||0p||c2(5«-i). 

Using estimates (12.231) , (I2.24p and Lemma 12. 3[ the existence of a solution to the fixed point 
problem (|2.22p can then be obtained by contraction mapping principle in the ball of radius 
2co ||0p|lc2(5«-i^ in the space Lg^g,{{R^ — {0}) ® SpanjxUoo}), provided that p is chosen larger 

than (but close enough to) Let us denote by Vp this fixed point. 

Since S > 1 — N, we have \vp\ << Up near and hence the solution u := (1 — x) Up + Vp is 
singular and positive near 0. We now prove that u :~ (1 — x) ""p + '^p is also positive at infinity. 
Indeed, the function Vp can be written as 

Vp = Vp + apX Uoo 
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where, according to formula (|2.2H) . ap can be computed as 

ap{N-l) 0Nda^ {A{1 - x) Up + \{l - x) Up + Vp\P) dx = \{1 ~ x) Up + Vp\P dx. > 

This implies that ap > 0, and by the maximum principle, it is now easy to check that u > in 
^+ ^ {0}- This completes the proof of Proposition ll.il □ 



2.3. Some open questions. The results of proof of Propositions 11.11 and 11.21 and their parallel 
with the radial case and p close to j^^, to which we recall ODE phase plane analysis applies, 
lead us naturally to several questions concerning existence of solutions of Am + = on the 
punctured half space M.^ — {0} with boundary data. We list some of them next. 

Question 1. We believe that the solution ui which has been obtained in Proposition 1 1 . II for p 
close to should actually exist for all p £ (^rjj ^^)' 

Question 2. When p = we believe that there exists a one parameter family of solutions of 

the form 

2 — JV 

u{x) = log \x\,6) 

where t i — > v{t, •) is periodic. This one-parameter family of solution corresponds to the well 
know periodic solutions for the singular Yamabe problem and also to Delaunay surfaces in the 
context of constant mean curvature surfaces. 

Question 3. When p > > 3, we believe that there exists a solution of Am + m^" = 

defined on M.^ which is identically equal to on dR^ and which is asymptotic to mq in (|1.10p 
at oo. This solution should correspond to the smooth radially symmetric solution of the same 
equation which is defined on the whole space and decays like at infinity, when p > ^^f. 

Question 4. Are there singular solutions when p > -^^j A^ > 4? In this regime separation of 
variables in general fails. 

Some partial answer to this question is given in [3]. 



3. The bounded domain case: proofs of Theorems 11.21 and 11.31 

The proof of our main results relies on two basic ingredients: one is the, already established, 
existence of the "basic cells" given by Propositions 11.11 and 11.21 which we will use to construct 
approximations to singular solutions. Another important ingredient, on which we elaborate in 
the next two subsections, is the analysis of invertibility of Laplace's operator, for right hand sides 
exhibiting a controlled singular behavior on a given embedded submanifold of dQ, in the same 
spirit to that of Lemma [2T3l Then we will use a fixed point scheme analogous to that in the proof 
of Proposition ll.il 

For notational convenience, we will assume in what follows that is actually a subset of R^, 
and that 5 is a smooth embedded submanifold of dil C with dimension k. We define 



N = n + k. 



16 



MANUEL DEL PINO, MONICA MUSSO, FRANK PACARD 



We start by setting up a suitable description of the space and Laplacian operator in natural 
coordinates associated to S. While the analysis below is done for A; > 1, it applies equally well 
to the point-singularity case k ^ 0, being actually simpler. 

3.1. Local coordinate system. In a neighborhood of a point p oi S let us choose coordinates 
yi, . . . ,yk on S. Next we choose sections Ei, . . . , En-i of the normal bundle of S in dil. We can 
define Fermi coordinates in some tubular neighborhood of S in dQ by using the exponential map, 

F{p, {xi,...,Xn-i))= Expf X, E,{p)) 

i 

for p in a neighborhood of p E S and (xi, . . . , x„_i) in some neighborhood of in K."^^. 

In these coordinates, it is well known that the induced metric gb on dQ can be expanded as 

gb = 5R"-i +9S + 0{\x\) 

where gs denotes the induced metric on S and x — (xi, . . . , Xn-i)- 

Finally, to parameterize a neighborhood of a point of dfl in Q, we denote by En the normal 
(inward pointing) vector field about and again use the exponential map to define 

G{q, Xn) = q + XnEn{q) 

for q in a, neighborhood of p in dfl and a;„ > in some neighborhood of 0. 

In these coordinates, it is well known that the Euclidan metric in Q can be expanded as 

5R„+fc = dxl +gb + 0{xn)- 

Collecting these two expansions, we conclude that in these coordinates the (Euclidean) Lapla- 
cian can be expanded as 

A = Ans + 0{\x\)V^+0{1)V (3.1) 

where x = (xi, . . . , x„) and Ans = ^R" +As is the Laplace Betrami operator on NS, the normal 
bundle of S in M". 

3.2. Analysis of the Laplacian in weighted spaces. We want to prove a result in the same 
spirit as that of Lemma in the current setting. To do this, we need to define weighted spaces 
on ri\S, which have a controlled blow up rate as S is approached. Unlike those in Lemma [2T3l we 
choose Holder spaces, which are more suitable to deal with linear perturbations which are second 
order operators. Let us define, for sufficiently small R > 0, half "balls" and "annuli" 

B+{R) ■.^{{p,x)eNS+ I \x\ e {0,R]} 

and 

A+{Ri,R2) ■.^{{p,x)eNS+ I \x\ e [Ri,R2]} 

B-i-{R) is roughly "half" of a tubular neighborhood of radius R of the manifold S, or just a ball 
in case that S reduces to a single points. We consider the following weighted space of functions 
defined on B+{R)\S. 



BOUNDARY SINGULARITIES FOR WEAK SOLUTIONS 17 

Definition 3.1. The space Cf°'{B+(R) \ S) is the space of functions u G Cf^"{B+{R) \ S)) for 
which the norm 

ll"llc^'°(S+(/?)\S) = sup r'^ \\u{-,r-)\\ce.^f^^+{r/2,r)) 

is finite. 

We consider now tiie problem 

Ansu = \x\-^f inB+iR)\S 

(3.2) 

u = OTLdB+{R)\S. 
We have the vahdity of the foUowing result. 

Lemma 3.1. Assume that (5 G (1 — n, 1). There exists a constant c > 0, independent of R> 0, 
such that, for each f G C^^'" {Bj^(R) \ S), there is a solution u — Gs,r (/) of problem (| j'.^p which 
defines a linear operator of f and satisfies the estimate 

M\cf''{B+iB.)\S) ^ ^\\f\\c'^-''{B+{B.)\S) ■ 



Proof. We only carry out the proof for i? = 1 since the general case follows by scaling. First we 
solve for each r G (0, 1/2) the problem 

Ansu = \x\-^f inA+(r,l) 

(3.3) 

u = on dA^{r, 1) 

and call Ur its unique solution. Maximum principle employed in a similar way as in Lemma |2.3[ 
taking into account expansion , yield the a priori bound 

\Ur\ < c||/||c«.°(B^(i)\s) 

where c = c(n, S) > 0. Then, elliptic estimates applied on geodesic balls of radius r centered at 
distance 2r from S give the following bound on the gradient of u 

iVu^l < c||/||co.»(B+(i)\s) 

for some c = c(n, S) > 0. Using Arzela's theorem, we conclude that, for a sequence of radii 
tending to 0, the sequence Ur converges to a function u which satisfies 

l«l < c||/|lc«-(i3+(l)\S) l^^l'^ 

and solves (|3.2p for R= 1. Again, elliptic estimates applied on geodesic balls of radius r centered 
at distance 2r from S yield the bound 

for some constant c = c{n, 6) > 0. Uniqueness of the limit u is easy to get and we leave it to the 
reader. The proof is concluded. □ 
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Next we will extend the previous result to the entire domain Q \ S. To do so, we consider a 
function 

j:n\S — > (0,oo) 

smooth, positive, which in the above defined local coordinates coincides with |a;| in a neighborhood 
of S in fj. This function will play the role of the function |a;| defined in B+{R) \ S. We define 
accordingly weighted Holder spaces as follows. 

Definition 3.2. We let the space Cg'"{fl \ S) be that of functions u e Cl'^ipl \ S) for which the 
norm 

\Mcl-''(n\S) = ll"llc*'°(B+(/?)\S) + ll"llc^-°(0\i3+(i?/2)) 

is finite. 

We consider now the problem 

Au = 7-2/ {nn\S 

(3.4) 

u = oiid^\S. 
We have the following result, extension of Lemma [3TT1 

Lemma 3.2. Assume that (5 G (1 — n^l). There exists a constant c > such that, for each 
f G C^'"(n \ S), there is a solution u — Gs (/) of problem \3.4\ which defines a linear operator 
of f and satisfies the estimate 

ll"llc,"'°(0\S) ^ C||/||c0.»(j^^5) . 

Proof. The proof follows from Lemma l3.11 expansion p.ip and a linear perturbation argument. 
First, we claim that the result of Lemma [3.11 remains true in B^{R) \ S* if the operator A^vs is 
replaced by A and if R is chosen small enough. Indeed, we have from p.ip and Proposition 13. II 

11/ - 7' (A - Ajvs) ° GsMMcl^iB+irnxs) ^ ci? ||/||cO,«(B^(^)\s)- 
The claim follows at once from a perturbation argument, provided that R is fixed small enough. 
We denote by Gs,r the right inverse for A in B+{R) \ S. 

We consider a cut-off function XR which is equal to 1 in _B+(i?/2) \ S and equal to in 
n\B+{R). We define 

/:=/-7'A(xflWi), 

where ui = GsMif)- Observe that this function is supported in \ i?+(i?/2). We have that 
/e C°'"(f7) and 

ll/llc«.°(fi) < 4f\\c°-''(n\r) 
for some constant c = c(n, S, R) > 0. 

Finally, we can solve 

Au2 — 7^^/ in Q 
M2 = on dil. 
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We have the bound 

Il'"2||c2.»(n) < c||/||co. = (f2\5). 

The desired result then foUows by letting the solution of (|3.4p he u = ui + U2- □ 

3.3. Proof of Theorems 11.21 and 11.31 We are now in a position to provide the proof of The- 
orem 11.21 and Theorem 11.31 The argument goes along the same lines as that in the proof of 
Proposition ll.il now with Lemma [5T^ plaving the role of Lemma [^751 

We recall that we are now assuming that is a domain in M^. We also write 

N = n + k. 

Proof of Theorem II. 2^ case p = . We assume that S is either a finite number of points of 
namely fc = 0, or an embedded A:-dimensional submanifold of dCl. For all e > small enough, 
we define 

Ue :=Xfle"~^"i(ek|,£a;„) 
where ui is the solution provided by Proposition 11.31 and Xft. is a cut-off function which equals 1 
in B+ (i?) \ S and in \ (2i?) . Here we fix i? > sufficiently small and use for x and Xn the 
meanings given in the previous subsections. In particular we have that = on dfl \ S. 

The problem we want to solve then reads 

A{u^ + v) + \u^ + = infl, 

V = oiidn\S, 

where we also require + v > in Q. Let us fix (5 G (1 — n, 2 — n]. By virtue of Lemma l3.2i we 
can rewrite this equation as the fixed point problem 

V = -Gs (l^i^Ue + +«|^)) . (3.5) 

We have the validity of the following fact: there is a constant cq = c(5, f2, S") > such that 

||72(A?/, + u|^)||cO.»(^^5) <co (log(l/£))T^, 

result that is a consequence of expansion p.ip and a direct computation using the asymptotic 
properties of ui in Proposition 1 1.21 

We restrict our attention to the case where 6 < 2 — n since 7^ (Awg -I- Ue^^ ) is bounded by a 
constant times near S, and S < 2 — n guarantees that this function belongs to Cg'°'{ri \ S). 

A second estimate wc can directly check is the following: Assume that d G {1 — n,2 — n] is 
fixed. There exists a constant c = c{6, ft, S) > such that 

\\l^{\u, + V2\^ - |u, +'yi|^)||co,»(a\s) < c(log(l/e))"^ 1^2 - «i|lc|-(o\s) 

for all V2,vi e Cg'°'{Cl \ S) satisfying 

ll^^llc,^'°(o\s) ^ 2co(log(l/e))V. 
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The above estimates allow an application of contraction mapping principle in the ball of radius 
2 Co (log(l/e))~2- in c|'"(f2\S') to predict existence of a solution to problem p.Sp . which we denote 
by Ve. 

Since S > 1 — n, we have \v^\ << near S and hence the solution u := u^ + v^ is singular along 
S and is positive near S. The maximum principle then implies that u > in Q. This completes 
the proof of Theorem 11.21 in the case p — ^^^j . □ 

The proof of Theorem II. 3^ case p = ^^zj- This proof uses similar arguments together 
with an induction process. By assumption, A is closed and contains a sequence of fc-dimensional 
submanifolds Si, i ^ N such that UiSi is dense in A. We define inductively the sequence of 
functions Ui which are solutions of 

Au, + uf^ = (3.6) 

in Q, satisfy Ui = on dfl — Uj^gSj and are singular along U*^o*5'j- Assume for example that 
Ui-i has already been constructed, then, we define 

Ui = +e""^Xri+i ui(eidist(-,S'i)) 

where ui is the solution provided by Proposition II. 3i is fixed small enough less than half the 
distance from Si to U^ZQS'j and > is small enough. Applying a perturbation argument as 

above, we can perturb Hi into a solution Ui = Ui + w,; of ()3.6p for some function Vi G C^'"(fi \ 
Wj^QSj). Taking Si small enough, we can ensure that 

||m»-w.-i||li(o) <2-' (3.7) 

||dist(-,91])2|u, <2-' (3.8) 

and 

\\J'v^\\L^m<2-' (3.9) 

where 7 = dist(-, S) and where 5 E (1 — n, 2 — n] is fixed. Clearly p.7p ensures that the sequence 
(iti)i converges in L^(il) to a function u. Moreover (|3.7p and (|3.8p imply that u is a weak solution 
of (jl.ip . Finally, (|3.9p implies that the nontangential limit of u at any point of S is equal to 
+00. □ 



Finally, we observe that, using Proposition II . II instead of Proposition II. 2) the results of Theo- 
rems [T2] and [T3] hold when p > is sufficiently close to -^^^ The only difference being that, 
in the proof of the result corresponding to the one of Theorem 11.31 in addition to the properties 
p.7p to p.9p which ensure the convergence of the sequence of solutions in the appropriate space, 
we may also ask that the sequence converges in W^''^{fl), for some q close enough to 1. The proofs 
are concluded. □ 
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